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Abstract 

The two-point function of the stress tensor operator of a quantum 
field in de Sitter spacetime is calculated for an arbitrary number of 
dimensions. We assume the field to be in the Bunch-Davies vacuum, 
and formulate our calculation in terms of de Sitter-invariant bitensors. 
Explicit results for free minimally coupled scalar fields with arbitrary 
mass are provided. We find long-range stress tensor correlations for 
sufficiently light fields (with mass m much smaller than the Hubble 
scale H), namely, the two-point function decays at large separations 
like an inverse power of the physical distance with an exponent pro- 
portional to m 2 /H 2 . In contrast, we show that for the massless case 
it decays at large separations like the fourth power of the physical 
distance. There is thus a discontinuity in the massless limit. As a 
byproduct of our work, we present a novel and simple geometric in- 
terpretation of de Sitter-invariant bitensors for pairs of points which 
cannot be connected by geodesies. 

1 Introduction 

Two key elements in our present understanding of cosmology [1, 2] are the 
existence of an early inflationary phase (which successfully explains the an- 
isotropics of the cosmic microwave background and the observed large-scale 
structure of the cosmos) and the present accelerated expansion of the uni- 
verse, which could be entirely driven by a non-vanishing cosmological con- 
stant according to present observations. The geometry of most inflationary 
models is close to de Sitter spacetime and so will be the late-time behavior of 
our expanding universe if it is currently driven by a cosmological constant. 



Therefore, the physics of de Sitter spacetime may be crucial to understand- 
ing the very early universe as well as its ultimate fate. 

The standard theoretical analysis relies on the linearized calculation of 
primordial cosmological perturbations generated in inflationary models [3], 
and the absence of large quantum corrections to the background dynam- 
ics of a spacetime driven by a cosmological constant. However, a growing 
number of recent studies [4, 5, 6, 7, 8] claim that nonlinear contributions 
from cosmological perturbations in inflationary models could actually be 
important. It has also been argued that higher-order radiative corrections 
from graviton loops could give rise to a secular screening of the cosmologi- 
cal constant for pure gravity in de Sitter spacetime [9, 10]. Such nonlinear 
effects are dominated by infrared modes (copiously generated by the ex- 
ponential expansion) and should be fairly insensitive to the short-distance 
behavior of a full theory of quantum gravity. However, when quantizing 
metric perturbations and considering their loop contributions, especial care 
is needed to make sure, by introducing appropriate diffeomorphism-invariant 
observables, that any nontrivial results obtained do not correspond to gauge 
artifacts [11, 12, 13, 14, 15, 16]. 

Our interest centers on the quantum metric fluctuations around exact 
de Sitter spacetime. At the linear level, these have been studied by several 
authors [17, 18, 19, 20]. The two-point function of the linearized Riemann 
tensor R a £ d (which is gauge-invariant around any maximally symmetric 
space) can be read from the results of Ref. [21], and it decays like the fourth 
power of the physical distance between the two points. Here we would like 
to go beyond the linear level and study the effects of matter loops on the 
quantum metric fluctuations. A key ingredient in this respect is the two- 
point function of the stress tensor operators of the matter fields. Calculating 
this object in de Sitter spacetime, and in a de Sitter-invariant vacuum, will 
be the main goal of this paper. This will be done using the formalism of 
de Sitter-invariant bitensors [22] . Besides offering a powerful computational 
technique, within this approach the results are directly obtained in a man- 
ifestly de Sitter-invariant form. This is particularly convenient, since one 
of the aims of this line of research is to gain insight on the possibility of a 
breaking of de Sitter invariance due to infrared quantum effects, as suggested 
by Tsamis and Woodard [9, 10] or Polyakov [23] (see also Refs. [24, 25] for 
related work concerning the latter proposal). 

Explicit results for free scalar fields in the Bunch-Davies vacuum [26, 27] 
are obtained. We assume minimal coupling to the curvature, but the mass m 
of the field is arbitrary. We find that, when m is small, the correlations have 
a long range: the stress tensor two-point function decays at large separations 
like an inverse power of the physical distance, with an exponent proportional, 
in the limit of small masses, to m 2 /H 2 , where H is the Hubble constant. 
Hence, the correlations can decay extremely slowly for sufficiently light fields. 
The strictly massless case is more subtle. As is well-known, in this case the 
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Bunch-Davies vacuum does not exist, and in fact there is no de Sitter- 
invariant vacuum at all [28]. Nevertheless, as explained in the appendix, 
one can consider a vacuum "arbitrarily close" to being de Sitter- invariant. 
We calculate the stress tensor two-point function of the massless field in this 
vacuum, and obtain a de Sitter-invariant result. Interestingly, we find that 
in this case there are no long-range stress tensor correlations, which instead 
decay at large separations like the fourth power of the physical distance. 
Such a discontinuity of the massless limit may seem surprising at first sight, 
but in fact it can be understood in simple terms. The key point is that in the 
massive case there is a contribution from the mass term in the stress tensor 
that does not vanish in the limit m — > 0, due to the infrared divergence of 
the Bunch-Davies vacuum. 

Finally, we would like to highlight a simple geometric interpretation of 
de Sitter-invariant bitensors for pairs of points that cannot be connected by 
geodesies which is presented in Sec. 4.1 and, to the best of our knowledge, 
has not been pointed out before. The basic de Sitter-invariant bitensors, 
in terms of which all other invariant bitensors can be constructed, are the 
geodesic distance between two points (i.e., the length of the shortest geodesic 
that joins them) and its first two covariant derivatives, which also have a 
simple geometric interpretation. However, in de Sitter spacetime there are 
pairs of points that cannot be connected by a geodesic. In Ref. [22] it 
was shown that, for these pairs of points, one can still define the geodesic 
distance and its covariant derivatives by analytic continuation, but no ge- 
ometric interpretation of the resulting objects was provided. We give this 
interpretation in Sec. 4.1. 

The paper is organized as follows. We start by reviewing in Sec. 2 the 
basic geometric properties of de Sitter spacetime that will be used in the rest 
of the paper (including a simple proof of the existence of pairs of points that 
cannot be connected by geodesies). The definition of bitensors and their 
general properties are introduced in Sec. 3 as well as those of maximally 
symmetric bitensors. The latter are studied in more detail for the particular 
case of de Sitter spacetime in Sec. 4, and a novel simple geometric inter- 
pretation for geodesically disconnected points is presented. These tools are 
employed in Sec. 5 to write the stress tensor two-point function for space- 
like separated points in de Sitter-invariant form, and explicit expressions for 
the case of a free minimally coupled scalar field are derived. Furthermore, 
closed results in terms of elementary functions are given in Sec. 6 for the 
particular cases of small masses, vanishing mass, and large separations. The 
extension to timelike separated points is provided in Sec. 7. We discuss the 
implications of our results for the quantum metric fluctuations around de 
Sitter spacetime in Sec. 8. The subtleties concerning the vacua of massless 
fields in the case of minimal coupling are discussed in the appendix. 

We use the (+,+,+) sign convention of Ref. [29] and the abstract index 
notation of Ref. [30]. Throughout the paper we work with natural units 
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Figure 1: Hyperboloid representing de Sitter spacetime. The point x and its 
antipode, a(x), are also represented, as well as a spacelike geodesic joining 
them. 

(h = c = 1) such that the Hubble constant equals one, unless explicitly 
stated otherwise. 



2 De Sitter spacetime 

We start by recalling the definition of de Sitter spacetime. Consider the 
(D + l)-dimensional Minkowski spacetime, A4 D+1 , with metric t\ab- The 
D-dimensional de Sitter spacetime, dSp, is the set of position vectors X A G 
Ai D+l that satisfy the equation 



VabX a X b 



H 



-2 



(1) 



where H is a parameter with dimensions of mass called the Hubble con- 
stant. This equation defines a hyperboloid, represented in Fig. 1. From 
this definition we see that dSp is the set of points of J\A D+1 that lie at a 
constant distance from the origin, and thus it is the Minkowskian analog of 
the sphere, S D . 

The isometry group of dSr> is formed by all the isometries of A4 D+l that 
leave the hyperboloid (1) invariant. These are the Lorentz transformations. 
Thus, the isometry group of dSp is the Lorentz group, 0(D,1). Since 
0(D, 1) is generated by D{D + l)/2 Killing vectors, which is the maximum 
number of Killing vectors that a manifold can have, dSr> is a maximally 
symmetric space, like S D . The group 0(D, 1) is also called the de Sitter 
group. 
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It can be shown that the metric of dSp, g a b, is a solution of Einstein's 
equations, 

G ab + ^9ab = 0, (2) 

where G a b is the Einstein tensor, and 

A = ' P - 1 » 2 ' P - 2 '^ (3) 

is the cosmological constant. Henceforth we will set H = 1. The dependence 
on H can be restored throughout by dimensional analysis. 



2.1 Points not connected by geodesies 

In dSo, some spacelike separated points cannot be connected by a geodesic. 
In order to see this, let us first introduce the biscalar 

Z(x,x') = VAB X A (x)X B (x'), (4) 

where x,x' £ dSp- Note that Z is invariant under the de Sitter group: 
if a is an isometry, i.e., X A {a{xj) = A A B X B (x), where is a Lorentz 
transformation, then Z(a(x),a(x')) = Z(x,x f ). It is also useful to introduce 
the Minkowskian distance, d(x,x'), defined by 

d 2 (x,x') = r, AB [X A (x) - X A (x')][X B (x) - X B (x% (5) 

That is, d(x, x') is the length of the straight line in the embedding Ai D+l 
that joins x and x' . Z and d are related through 

d 2 = 2(l-Z), (6) 

as follows from Eq. (1). Therefore, Z = 1 when x' is on the light cone 
of x, Z > 1 when the points are timelike separated, and Z < 1 when the 
points are spacelike separated. It should also be noted that, as the squared 
Minkowskian distance between two points on the hyperboloid can take any 
value, Z £ (— oo, oo). 

Consider now the spacelike curve represented in Fig. 1, passing through 
the point x. This curve is a geodesic, because it is a straight line on a 
cylinder tangent to the hyperboloid. The Minkowskian distance between 
x and any other point x' in this geodesic satisfies < d(x,x') < 2. This 
implies 

-1 < Z{x,x) < 1. (7) 

Now, any other spacelike geodesic passing through x can be obtained from 
this one by an isometry that leaves x invariant. This is analogous to the 
case of the sphere, where any geodesic (meridian) passing through the north 
pole is obtained from a given one by a rotation that leaves the north pole 
invariant. This means that any point y' connected to x by a spacelike 
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Figure 2: The conformal diagram of de Sitter spacetime. The points in the 
shaded regions are geodesically connected with x. 

geodesic satisfies y' = a(x') for some x' in the geodesic represented in Fig. 1, 
and where a is an isometry such that cr(x) = x. De Sitter invariance of Z 
then implies Z(x,y') = Z(a(x), a(x')) = Z(x,x'). Consequently, equation 
(7) holds for any point x' connected with x by a spacelike geodesic. In other 
words, spacelike separated points x, x' such that Z(x,x') < —1 cannot be 
joined by a geodesic. 

There is no similar restriction for timelike geodesies, as there are the- 
orems that guarantee the existence of geodesies between any two timelike 
separated points when the spacetime is globally hyperbolic [30]. The re- 
gions of dSp that are geodesically connected to a point x are depicted in 
the conformal diagram of dSp in Fig. 2. 

2.2 The antipodal transformation 

We will devote special attention to one particular isometry of dSu called 
antipodal transformation, that we denote by a. This operation sends a 
point x to its antipodal point a(x), defined by X A {a{x)) = — X (x). An 
example is given in Fig. 1. This isometry has some special features, namely 
that it commutes with the rest of elements of the de Sitter group, o{a{x)) = 
a(a(x)), and that its square is the identity, a{a{x)) = x. On the other 
hand, it follows directly from the definition of the antipodal transformation 
that 

Z(a(x), x) = — Z(x, x). (8) 

Consequently, when x' is on the light cone of a(x), then Z(x, x') = —1. This 
is shown in Fig. 2. 
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3 Bitensors 



In this paper we want to compute the stress tensor two-point function of a 
quantum field. What kind of object is this, from a geometric point of view? 
It is not a tensor, because it takes values on every pair of points of the 
spacetime, and it has indices associated with each of the points. It is rather 
a bitensor. Bitensors were first introduced by Synge [31]. In this section we 
will give a precise definition of what bitensor means, and summarize some 
results that will be important for us. 

Let M be a manifold. Consider two points x, x' £ M, and let V x and V x i 
be the tangent spaces at x and x' respectively. A bitensor, T a ...^, of type 
(0, k)(0, k') over V x and V x > is a multilinear application 1 

T a ...v '■ Vx x • • • x V x x V x > x • • • x V x i — ► C 

That is, given k vectors at the point x and k' vectors at the point x' , T a j,i 
produces a complex number. As we see, the definition of a bitensor resembles 
that of an ordinary tensor. In particular, a real-valued bitensor of type 
(0, k) (0, 0) is a tensor of type (0, k) . 

A bitensor field, T a . b', is a rule that assigns a bitensor over V x and V x / , 
T a ...b'{ x -, %'), to each pair of points (x, x') £ M x M. Note that a real- valued 
bitensor field of type (0, fc)(0, 0) is not a tensor field, but a rule that assigns 
a tensor over V x to each pair of points (x, x'). In spite of this, we will refer 
to bitensor fields simply as bitensors, as is commonly done with ordinary 
tensor fields. 

A bitensor can be differentiated with respect to either of its arguments. 
The covariant derivative of a bitensor, T a ..jy , of type (0, k) (0,k') with respect 
to x, which is denoted by V c T a ...&', is obtained by (i) fixing the value of x'\ 
and (ii) taking the covariant derivative of the resulting object as if it were a 
tensor of type (0, k). The covariant derivative with respect to x' is defined 
analogously. 

Bitensors of type (0,0)(0,0) and (0, 1)(0, 1) are also called biscalars and 
bivectors respectively. Bitensors are central mathematical objects in quan- 
tum field theory, because any two-point correlation function is a bitensor. 
For example, the propagators of the scalar, the photon and the graviton 
fields are a biscalar, a bivector and a bitensor of type (0,2) (0,2) respec- 
tively. 

3.1 Invariance of a bitensor 

The notion of invariance of a bitensor is, again, closely related to the well- 
known notion of invariance of a tensor. Consider a diffeomorphism a : M — ► 

1 One can define more general bitensors, of type (j, k)(j' , k'), which have superindices as 
well as subindices. Here we restrict ourselves to bitensors of type (0, k)(0, k') for simplicity. 
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M, that assigns a point a(x) to each point x, and let a* be the pushforward 
associated with this diffeomorphism. That is, if the vector v a £ V x is tangent 
to the curve x(A) at the point x, then the vector a*v a £ is tangent to 
the curve a(x(\)) at the point a(x). A bitensor T a ..v is said to be invariant 
under a if 

T a ... b ,(a(x),a(x'))a*v a . . . a*w b ' = T a ... 6 ,(x, x')v a ...w b ' (9) 

for all pairs of points (x, x'), and for any v a , . . . £ V x and . . . , w b £ V x i. An 
important property concerning invariance of bitensors is that, in the case 
where a is an isometry, the covariant derivative of a bitensor invariant under 
a is also invariant under a. 

3.2 Maximally symmetric bitensors 

Consider the case where M is maximally symmetric. A bitensor is said 
to be maximally symmetric if it is invariant under all the isometries of M. 
The following bitensors are maximally symmetric [22]: the distance, fi(x, x'), 
along the shortest geodesic joining x and x', also called the geodesic distance; 
the unit vectors, n a (x, x') and n a '(x, x'), tangent to the geodesic at the points 
x and x' respectively, pointing outward from it; the parallel propagator, 
g a b'( x , x '), which parallel-transports a vector from x' to x along the geodesic; 
and the metric tensors, g a b(x,x') and g a i(,'(x,x'), at the points x and x' 
respectively. Furthermore, it was shown in Ref. [22] that these bitensors form 
a complete set, in the following sense: any maximally symmetric bitensor is a 
linear combination of products of n a , n a r,g a b, g a iy and g a b', with coefficients 
that depend only on fx. 

It is worth being more precise with the definition of fi,n a and n a i. If 
we denote the shortest geodesic joining x and x' by x(X) and its tangent 
vector by v a (X), and we take x(0) = x and x(l) = x', the geodesic distance 
between x and x' is 

ti(x,x')= / 1 dA[ ffa6 (x(A))^(A)^(A)] 1 /2. (io) 
J o 

Note that /j,(x,x') is a real number for spacelike separated points, and an 
imaginary number for timelike separated points. The sign prescription for 
the square root is chosen in such a way that fi(x, x') > for spacelike 
separated points, and fi(x,x') = i£(x,x'), with £(x,x') > 0, for timelike 
separated points. The bitensors n a and n a i are defined by the equations 

n a = V a /X n a , = V a 'fi. (11) 

The geometric interpretation of n a (x, x') and n a '(x, x') is obtained by differ- 
entiating Eq. (10). One finds that they are the unit vectors tangent to the 
geodesic at the points x and x', pointing outward from it, and multiplied by 
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the imaginary unit in the case where x and x' are timelike separated. This 
implies n a n a = n a m a = 1 for both spacelike and timelike separated points. 

Some useful properties of n a , n a i and g ab >, which follow directly from the 
definition of these bitensors, are 

g ab /(x,x')n b ' \x,x') = -n a (x,x') 

9a ( x : x)g db {x\ x) = g a b{x). (12) 

From the last of these equations, and the fact that parallel transport leaves 
the scalar product invariant, one can prove the additional relation 

g ab >(x,x f ) = g Va (x',x). (13) 

The theorem of Ref. [22] can be applied to the computation of the covariant 
derivatives of n a , n a > and g ab i, which are also maximally symmetric bitensors. 
The corresponding jU-dependent coefficients were also computed in Ref. [22] 
for all maximally symmetric spaces. In the particular case of dSo, the result 
is 

V a n b = cot n (gab - n a n b ) 
V a n b / = - csc/x (g aV + n a n b <) 

^a9bd = (csc^ - cot //) (g ab n c/ + g ac >n b ). (14) 

4 De Sitter-invariant bitensors 

The bitensors /j,,n a ,n a > and g ab > of dSn are defined in the domain 

Q = {(x, x') G dSo x dSo such that 3 a geodesic joining x and x'}. 

As we know from Sec. 2, this does not cover the whole dSo x dSp: only 
the pairs of points with Z > — 1 belong to Q, where Z has been defined in 
Eq. (4). However, as was shown in Ref. [22], the definition of fi,n a ,n a ' and 
g ab i can be extended to pairs of points outside Q in such a way that all their 
properties are preserved. Let us review the argument of Ref. [22]. 

For spacelike geodesies, — 1 < Z < 1, the geodesic distance /x of dSo is 
given by the equation 

H = cos" 1 Z, (15) 

where cos -1 is the principal value of the inverse cosine, which is single- 
valued. For example, cos -1 1 = 0, and cos _1 (— 1) = n. For the geodesic of 
Fig. 1, this equation is easily seen from Fig. 3. Other spacelike geodesies are 
obtained from this one by an isometry, and therefore, by de Sitter invariance 
of n and Z, Eq. (15) is true for all spacelike geodesies. 

Using a similar argument, it can be shown that a similar equation holds 
for timelike separated points, Z > 1: 

fj, = cos^ 1 (Z — ie), (16) 
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Figure 3: The geodesic of Fig. 1, seen from above. The geodesic distance 
fi(x, x') between x and x' is the angle between the position vectors X A (x) 
and X B (x r ). 

where e is a positive real infinitesimal. This —ie term is added because the 
function cos -1 Z has a branch cut at Z < — 1 and Z > 1 along the real axis, 
across which its imaginary part changes sign. 

Eqs. (15) and (16) serve to define \x outside Q by analytic continuation. 
Let us label the points of dSn with their corresponding position vectors in 
the embedding A4 D+1 . The geodesic distance, /j,(X , Y A ), is only defined for 
pairs of points (X A ,Y A ) such that Z(X A ,Y A ) > -1. However, Z(X A ,Y A ) 
is an analytic (polynomial) function of (X A ,Y A ) for all complex values of 
X A and Y A . Furthermore, cos -1 Z is an analytic function of Z in the whole 
complex plane except the branch cut. Therefore, cos -1 Z is the analytic 
continuation of to complex values of X A and Y A . We can define fi for 
pairs of points with Z < — 1 as 

fi = cos _1 (Z + k), (17) 

where the sign of the ie term has been taken positive for convenience. One 
could also adopt the opposite choice, in which case the imaginary part of 
the geodesic distance \i outside Q would have the opposite sign. 

The definition of n a and n a > is extended outside Q by substituting (17) 
in (11). The definition of the parallel propagator can also be extended using 
the second equation in (14), which gives g a y in terms of fi, n a and n a i. Since 
these definitions have been obtained by analytic continuation, all the results 
concerning fj,, n a , n a i and g a y that we have presented in the previous section 
hold everywhere in dSp x dSo- 

4.1 Geometric interpretation 

What was not emphasized in Ref. [22] is that the bitensors fi,n a ,n a i and 
g a b> also have a geometric interpretation outside Q. We will explain it here. 
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Figure 4: Example of unit vectors n a and n a i for points x and x' that are 
not connected by a geodesic. 

The principal value of the inverse cosine has the property 

cos _1 Z = 7r-cos _1 (-Z). (18) 

On the other hand, as we know from Sec. 2, Z(a(x), x') = —Z(x, x'). There- 
fore, the equation above and (15)-(17) imply 

fj,(x, x') = 7T — /j,(a(x), x') (19) 

throughout dSo x dSo- For pairs of points in the geodesic of Fig. 1, this can 
also be seen from Fig. 3. Now, for (x,x') ^ Q we have Z(x,x') < —1, which 
implies Z((a(x),x') > 1. Therefore, a(x) and x' are geodesically connected 
and timelike separated. This can also be seen in Fig. 2. This means that 
the geodesic distance for geodesically disconnected points is ir minus the 
length of the shortest geodesic joining a{x) and x' . Since this geodesic is 
timelike, the geodesic distance is a complex number outside Q: \i = ir — i£, 
with 

Substituting (19) in (11), we obtain 

n a (x,x') = -n a (x,a(x')) n a ,{x,x) = -n a '(a(x),x'), (20) 

where we have used that fj,(a(x),x') = fx(x,a(x')). This is because \i is de 
Sitter-invariant, and the antipodal transformation is an isometry such that 
a(a(x)) = x. Eq. (20) means that, for (x,x') ^ G, n a (x,x') is the unit 
vector tangent at x to the geodesic joining x and a(x'), and n a i (x, x') is the 
unit vector tangent at x' to the geodesic joining a(x) and x', both of them 
pointing inward. This is shown in Fig. 4. 

Finally, after substitution of (19) and (20) in the second equality of (14), 
we obtain for the parallel propagator 

9ab' = -9ab> ~ 2n a n a <, (21) 
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where the bivector g a y is defined by the equation 

g~ab'(x, x)w b ' = g a b'(x, a(x'))a*w b ' (22) 

for all vectors w b ' £ V x >. In order to obtain Eq. (21), we have used the 
property ny(a(x), x')w b = ny(x, a(x'))a*w b , which holds for all vectors 
w b G V x i , and which is a direct consequence of de Sitter invariance of ny and 
the properties of the antipodal transformation. For (x,x') ^ Q, g a y(x,x') 
is the operator that "pushes forward" a vector from x' to a(x'), and then 
parallel transports it from a(x') to x along the geodesic joining these points. 

4.2 A comment on the Minkowskian distance 

In the previous subsection we have found the geometric interpretation of 
the bitensors n,n a ,n a > and g a y for pairs of points that are not connected 
by a geodesic. In the particular case of the geodesic distance [/,, note that, 
although the geometric interpretation is simple, fi is not the length of any 
spacelike curve (remember that x and x' are spacelike separated when they 
are not geodesically connected). Therefore, it cannot be taken as a notion 
of distance. In order to interpret physically our results below for spacelike 
separated points, it will be more convenient to express them in terms of the 
Minkowskian distance, defined in Eq. (5). The reason is that, for x and 
x' spacelike separated, this biscalar coincides with the so-called "physical 
distance" or "proper distance" associated with the flat slicing of dSo- Let 
us see this. If we label the points in dSo with the coordinates {t, x 1 }, with 
i = 1, . . . ,D — 1, defined by 

X° = - sinht - ^5ijX l x j e t 
X i = xV 

X D = cosht - ifcjzVe*, (23) 

then the line element takes the form 

ds 2 = r) AB dX A dX B = -dt 2 + e 2t 5 ij dx i dx j , (24) 

so that the surfaces of constant t are flat. Moreover, any two space- 
like separated points can be brought to the same flat hypersurface by an 
isometry. Now, for x,x' £ E$ we have 

d 2 (x,x') = ePSijtfix) - - x j (x% (25) 

as can be straightforwardly seen by using Eqs. (5) and (23). Hence, for 
x, x' G £t the Minkowskian distance coincides with the length of the straight 
line of Et that joins x and x'. This is what in cosmology is called "physical 
distance" or "proper distance" . 
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5 Stress tensor fluctuations in de Sitter spacetime 



Consider a quantum field in dSo, whose stress tensor operator is T ab , in a 
vacuum state |0) that we assume de Sitter-invariant. We want to compute 
the two-point function of the stress tensor operator, 

F abdd ,{x,x') = (0\T ab (x)T c , d ,(x')\0) - (0|T a6 (x)|0)(0|T c ^(x')|0), (26) 

which we will also refer to as the stress tensor fluctuations. For the moment 
we will assume that the points x and x' are spacelike separated. What we 
have seen so far tells us a lot about this two-point function before doing 
any computation, even before specifying the quantum field theory under 
consideration. Indeed, if we choose |0) to be a de Sitter- invariant state, 
the expectation value of any observable must be de Sitter-invariant. If the 
points x and x' are spacelike separated, causality requires 

[T ab (x),T c , d ,(x')}=0, (27) 

which implies that F abc / d i is the expectation value of a Hermitian operator, 
i.e., an observable (the product of two Hermitian operators that commute is 
Hermitian). Therefore, when x and x' are spacelike separated, F abc i d i is a de 
Sitter-invariant bitensor. According to the theorem stated in the previous 
section, it must be a linear combination of products of n a , n a i ,g a b,9a'b' an d 
g ab >, with coefficients that depend only on fi. On the other hand, for spacelike 
separated points Eq. (27) implies 

Fabc'd'(x,x') = F c , d , ab (x',x). 

And, since the stress tensor is symmetric, we have 

Fabc'd' = F bac / d ' = F abd ' c '- 

The most general bitensor of type (0, 2)(0, 2) (i.e., with the index structure 
of F abc 'd') that is a linear combination of products of n a ,n a > , g a b, 9a'b' and 
g ab >, with coefficients that depend only on fi, and that satisfies the above 
constraints is 

Fabc'd' = P{v)n a n b n c m d ' +Q{li){n a n b g c i d i +n c ,n d /g ab ) 

+R(li)(n a n c ,g hd , + n b n d ,g ac i + n a n d ,g bc , + n b n c ,g ad ,) 
+S(fi)(g ac 'gbd' + gbc'gad') + T(/J,)g ab g c , d >. (28) 

Moreover, the stress tensor conservation law implies 'V a F a b c 'd' = 0. Imposing 
this equation on (28), and using (14), we obtain 

P' + Q'- 2R! + {D- l)P cot n - 2Q esc n 

-2R[(D - 2) cot n - D esc fi\ = 

Q' + T / + (D-l)Qcot/i-2i?csc/i-2S'(cot/i-csc^) = 

R' - S' - Qcscn + DRcotfi- DS(cot [i- cscfi) = 0, (29) 
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where the prime stands for derivative with respect to fx. Eqs. (28) and 
(29) are valid for any quantum field theory in dSp, assuming that |0) is de 
Sitter-invariant and that the points are spacelike separated. 

There is yet another constraint on the functions P, Q, R, S and T, which 
comes from requiring the stress tensor two-point function to be unambiguous 
for antipodal points, x' = a(x). The unit vector n a is not well-defined when 
x' = a(x), because there is an infinite number of geodesies joining x and 
a(x), all of them with the same length, n(x, a(x)) = ir. This is analogous to 
what happens in the sphere: the north and the south poles are joined by an 
infinite number of meridians, all of them with the same length. The parallel 
propagator is also not well-defined, but it can be split into an ambiguous 
and an unambiguous part through Eq. (21). Using this equation, (28) can 
be rewritten as 

Fabdd' = (P + 8S - 8R)n a n b n c m d ' + Q{n a n h g c i d > + n c m d ,g ab ) 

+(2S - R){n a n c ,g bd > + n b n dl g ac! + n a n d ig bc , + n h n d g ad i) 
+S {g ac '9bd' + 9bc'9ad') + Tg ab g c/d > . (30) 

The bitensors appearing in this equation are all well-defined for antipodal 
points except for n a . In order for the stress tensor two-point function to be 
unambiguous when x' = a(x), the coefficients of the terms containing n a in 
the equation above should vanish at fx = w, 

P(tt) + 8S(vr) - 8-R(ir) = Q(tt) = 2S(tt) - R(tt) = 0. (31) 

This equation is also valid for any quantum field theory in dSp, assuming 
that |0) is de Sitter-invariant. 

5.1 Stress tensor fluctuations of a free scalar field 

Let us be more concrete and consider a particular quantum field theory in 
dSp, namely that of a free, minimally coupled scalar field (f>, with mass m. 
The corresponding equation of motion is the Klein-Gordon equation 

(V a V a - m 2 )(t> = 0. (32) 

The stress tensor operator of this field can be written using point splitting 
regularization as 

T ab (x)= lim V ab (x,x")(j)(x)^(x"), (33) 

with 

"Dab = 5 b c "V a V c » - -5a65 cd "V c V d » - -m 2 g ab . (34) 

In order to check that this expression coincides with the standard one [27], 
one has to take into account that the parallel propagator becomes the metric 
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in the coincidence limit, i.e., g a w>(x,x") — > g a b(%) when x" — > x. The stress 
tensor two-point function of the free scalar field is obtained by substituting 
(33) in (26), 

F abc , d ,(x,x') = lim V ab {x,x")V dd ,{x',x'") [(0|<K^(z>(zO<Kz"')|0) 

x"— >x 
x"'^x' 

-(O|0(x)0(x")|O)(O|0(x')0(x'")|O)]. (35) 

Since the field <p is free, and |0) is a vacuum state, the four-point function 
appearing in this equation can be written in terms of the Wightman function, 
G + (x,x') = (O|0(x)0(xO|O), 

(O|0(x)0(x")0(x , )0(x" , )|O) = G + (x,x")G + (x',x"') + G + (x, x')G + (x", x'") 

+G + (x,x"')G + (x",x'). (36) 

This equation is obtained by expanding the field (j) in terms of creation and 
annihilation operators. Note that, in the particular case where the product 
of fields on the left-hand side is time-ordered, this equation also follows from 
Wick's theorem. Substituting (36) in (35), and taking the limit x" — > x and 
x'" x', we find 

Fabc'd' = V a V c /G + V 6 V d /G+ + V a V d /G + V 6 V c /G+ 

-9ab {V e V c ,G + V e V d ,G + + m 2 V d G + V dl G + ) 

-9c'd> (v a V e ,G+V 6 V e 'G + + m 2 V a G+V 6 G + ) 
V e V P /G + V e V e 'G+ 



1 

9ab9c'd> 



2 



+m 2 (v e G + V e G + + V e /G+V e 'G+) + m 4 G+ 2 ] . (37) 



Note that, after substitution of (36) in (35), there has been a cancellation of 
terms that diverge in the limit x" — > x, x'" — ► x', so that this result is finite 
as long as x is not on the light cone of x'. 

At this point we impose that the vacuum |0) is de Sitter-invariant, and 
that the points x and x' are spacelike separated. For spacelike separated 
points, the Wightman function is the expectation value of an observable, so 
that it must be a de Sitter-invariant biscalar. Thus it must have the form 

G + = G(n). (38) 

On the other hand, the Wightman function satisfies the Klein-Gordon equa- 
tion (32). Using (38) and (14), this equation translates into 

G"(n) + (D - 1) cot nG'(n) - m 2 G{fi) = 0. (39) 

Imposing as boundary conditions that G(fi) diverges only at the coincidence 
limit, fx — > 0, and that it does so in the same way as the standard Wightman 
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function in Minkowski spacetime, the only solution is [32] 



G(fi) = c nhD F(h + ,h_; ^; i±^), (40) 

where F denotes the hypergeometric function, h± and c„ h D are constants, 

1 



(D - 1) ± ^(D — l) 2 — 4m 2 
(47r) D / 2 r( 



C m ,D = ~ Nn/owg , , (41) 



and we are using that Z = cos /x, as follows from (15) and (17). One can see 
[27] that there is a vacuum whose Wightman function for spacelike separated 
points is (40). It is the so-called Bunch-Davies vacuum. In what follows 
we will be interested in the stress tensor fluctuations in the Bunch-Davies 
vacuum. 

Substituting (38) in (37), and using (14), we find that the stress tensor 
two-point function in the Bunch-Davies vacuum has the form (28), as corre- 
sponds to a de Sitter-invariant state. The associated /^-dependent functions 
are 

P = 2G\ 

Q = -G 2 1 + 2G 1 G 2 -m 2 G' 2 
R = G1G2 
S = G\ 

T = 1 2 G 2 1 -G 1 G 2 + D -^G 2 2 + m 2 G' 2 + ^m 4 G 2 , (42) 
where, again, the prime stands for derivative with respect to fx, and 

G^n) = G"(n)-G'(n)cacn 

G 2 (fx) = -G'(/i)csc/A- (43) 

The functions (42) satisfy the stress tensor conservation equations (29), as 
can be checked by using the Klein-Gordon equation (39). They also satisfy 
the conditions (31), as follows from basic properties of the hypergeometric 
function. 



6 Particular cases 

Eqs. (42) give the stress tensor two-point function of the free, minimally cou- 
pled scalar field in terms of the Wightman function (40) and its derivatives 
with respect to fx. Of course, this result is not very explicit. In this section 
we will derive explicit expressions for the stress tensor two-point function 



16 



in the following particular cases: small masses, strictly zero mass, and large 
negative values of Z. 

The results will be given in terms of the biscalar Z. In order to interpret 
them, recall that Z is related to the Minkowskian distance d through Eq. (6). 
Remember also, from Sec. 4.2, that the Minkowskian distance coincides with 
the physical distance when the points are spacelike separated. 



6.1 Small masses 

Let us assume that the mass of the field is small, m <C 1. Remember 
that we are using units in which H = 1. In a general system of units, our 
assumption is that the mass is much smaller than the Hubble constant. This 
is the situation one encounters in inflationary models. 

We begin by studying the behavior of the Wightman function (40) in the 
limit of small masses. Using the definition of the hypergeometric function 
[33], Eq. (40) can be rewritten as follows for |1 + Z\/2 < 1: 

c(u)- 1 y r(ft + + n)r(ft. + n) i (i + zy 

For other values of Z, the Wightman function is obtained by analytic con- 
tinuation. The dependence on m is implicit in the parameters h±, Eq. (41). 
Expanding the gamma functions in the equation above in powers of m we 
obtain at leading order 

The Wightman function diverges in the limit m — > 0. This is the well-known 
infrared divergence of the Bunch-Davies vacuum. Since the divergent term 
is independent of fi, the derivative of the Wightman function G'(/x) will be 
finite in this limit. Indeed, expanding the derivative of (44) in powers of the 
mass, we obtain after some algebraic manipulations 

fc=0 ^ ' 

for an even number of dimensions D, where 

_ T(D/2) ( D-l \ 
Ck -J^)m\D/2 + k)- (47) 

Note that (46) is a sum over a finite number of terms, not an infinite series. 
Eqs. (45) and (46) have been derived assuming |1 + Z\/2 < 1. However, 
since the extension to other values of Z is obtained by analytic continuation, 
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these equations are actually valid for any Z < 1. This result generalizes to 
an arbitrary, even number of dimensions a previously obtained expansion of 
the Wightman function in powers of the mass in dS^ [34] , and reduces to it 
when we set D = 4. 

Now that we have written the Wightman function and its derivative 
with respect to ji in terms of simple functions, we can obtain an explicit 
expression for the stress tensor fluctuations. Substituting (45) and (46) in 
(42), we find that the functions P, Q, R, S have the form 



1 



(i-zy 



eV^" 



n=0 



0(m 2 ), 



(48) 



where we are denoting P, Q, R, S generically by F. For the remaining func- 
tion we obtain 



D-2 



Tin) 



{i-zy 



5> 



n=0 



l + Z 



+ 



1 



T(D) 



(4vr)^r(f) 



+ 0{m 2 ). (49) 



In these equations, the coefficients F n and T n are 



A? 



n + I 



Qn — 
Rn = 
S n = 



2 

n + l 
2 

n + I 



i=-p 

i 

- 4 E 

z=- 

^E 

l=—N 
N 

EC n+l C n-i 

l=-N 
J=-JV 



+ 1 



n — I 



+ 1 



2 

n — I 



+ 1 ) Cn+l Cn-l 
2 2 



c n—l 

2 2 



+ 1 )C n+l c ■n 



+ 1 



n — I 



+ 



D-4 



C n + l C n — l . 
2 2 



(50) 



where iV = min(n, D — 2 — n). 

The stress tensor two-point function of a free scalar field with small mass, 
for spacelike separated points, is given by Eqs. (48) and (49). It is finite in 
the limit m — > 0, unlike the Wightman function. This is because Eqs. (42) 
involve only derivatives of G, which are finite, with the only exception of 
T, which includes a term proportional to m 4 G 2 . The factor m 4 cancels the 
divergence of the Wightman function. 

An interesting aspect of the obtained stress tensor two-point function is 
its behavior at long Minkowskian distances. From Eq. (6), this corresponds 
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to large negative values of Z. Setting Z <C — 1 in (48) and (49), we obtain 



P, Q, R, S 



T 




(51) 



From Eq. (6), we see that the leading terms in P,Q,R,S decay like d~ 4 . 
More interestingly, the leading term in T does not decay to zero. In other 
words, the stress tensor two-point function in the limit m — > does not 
vanish at long distances. 

The origin of this behavior is, again, the term proportional to m 4 G 2 
in Eqs. (42). Indeed, the terms with derivatives of G are either 0(m 2 ), 
and thus subleading, or vanishing in the limit Z — > — oo, whereas the term 
proportional to m 4 G 2 is a constant of order 0(m°). 

It is worth comparing this result with the stress tensor fluctuations in 
Minkowski spacetime, Ai D . Note first that, since the theorem stated in 
Sec. 3 is valid for all maximally symmetric spaces, the stress tensor fluc- 
tuations in A4 D also have the form (28). The corresponding functions 
P, Q, R, S, T can be obtained by dimensional analysis. Indeed, these func- 
tions must have dimensions of (energy density) 2 , that is, (mass) 215 . But 
the theory of a free scalar field in M D , and in the limit m — > 0, has no 
dimensionful constants. Thus, by dimensional analysis the only possibility 
is 



Hence, the stress tensor fluctuations in the limit m — > do decay with the 
distance in M D . Finally, it is worth noting that, for small distances (Z — > 1) 
the functions (48) and (49) have the same form as in M D , Eq. (52), as it 
should be. 

6.2 Massless field 

We now consider a slightly different case, namely that of a strictly massless 
field, m = 0. As is well-known [28], in this case there exists no de Sitter- 
invariant vacuum. The vacuum states we choose for the massless field are 
the 0{D) -invariant vacua, |0) a , introduced by Allen and Folacci [34]. These 
states form a one-parameter family with a £ (0,oo). We will be interested 
in the limit a — > 0, because, although |0) Q are not de Sitter-invariant, there 
is a class of observables whose expectation values are de Sitter-invariant in 
this limit. This is explained in the appendix. 

The stress tensor two-point function of a free scalar field is given by 
Eq. (37), for any mass m and for any vacuum state. In the massless case 



P, Q, R, S, T oc n 



-2D 



(52) 
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this equation reduces to 



Fabc'd* = V a Vc>G + V b V d ,G + + V a V d ,G + V b V c ,G + 

-g ab VeVc>G + V e V d ,G + - g c , d ,VaVeiG + V b V e 'G + 

+ ^a^ c 'd'V e V e ,G+V e V e 'G+. (53) 

Note that all terms in this equation are products of V a V&'G + . Let us denote 
by F°L , d , the stress tensor two-point function in the vacuum |0) a . We want 
to compute it in the limit a — ► 0. Now, in the appendix it is shown that 

lim V a V 6 /G+ = lim V a V 6 /G+ D , (54) 

a— >0 m^O 

where G+ is the Wightman function of the massless field in the vacuum \0) a , 
and G~^ D is the Wightman function of the massive field in the Bunch-Davies 
vacuum. Therefore, the stress tensor two-point function of a massless field in 
the limit a — > can be computed as follows: (i) compute (53) for a massive 
field in the Bunch-Davies vacuum; and (ii) let m — ► 0. It is important to 
notice that this is not exactly the same as computing the stress tensor two- 
point function of a massive field in the limit m — > 0, which is what we have 
done in the previous subsection. 

The step (i) corresponds to substitute the Wightman function (38) in 
(53). The result has the form (28), with 

P = 2G\ 

Q = -G? + 2GiG 2 

R = G1G2 

S = G\ 

T = ^Gl-G^ + ^-^Gl (55) 

and where G± and G2 are given by (43). This equation is exactly the same 
as (42), but without the mass terms. We can now let m — > 0. Following the 
arguments of the previous subsection, we obtain 

where, as before, F denotes generically P, Q, R and S, and F n and T n are 
given by Eqs. (50). We have thus obtained for the massless field exactly the 
same result as in the limit m — ► 0, except for the constant piece of 
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in Eq. (49). Remember that such a constant piece arises because of a mass 
term, which is lacking here. 

Due to the absence of the constant piece, the behavior at large distances 
changes dramatically with respect to the case of small masses: at Z <C — 1, 
we have 

P,Q,R,S,T ~ Z- 2 . (57) 

Therefore, the stress tensor two-point function of the massless field vanishes 
at long distances. There is a "discontinuity" at m = 0: taking the limit 
m — > is not the same as setting m = 0. 

One might think that this discontinuity is due to our choice of states for 
both the massive and the massless fields, and that by choosing appropriate 
states the discontinuity would disappear. However, this does not seem to 
be the case. In Ref. [35], the expectation value a (0|T a b|0) Q was computed 
for all values of a. In the limit a — > the result is de Sitter- invariant and it 
disagrees with the limit m — > 0, similarly to what we have just seen for the 
stress tensor two-point function. Furthermore, at late times one has 

a(0|T a6 |0) a lim a (0|T ab |0) Q . 

That is, lim a _>o a(0|T a 5|0) Q is an attractor at late times. As is well-known, 
the same behavior occurs for the massive field in the Bunch-Davies vacuum. 
This means that, in the case of the expectation value of T &, the discontinuity 
at m = appears at late times regardless of the state of the quantum field. 
It is plausible that something similar happens for the stress tensor two-point 
function. 

It is important to note that we are comparing the massless, minimally 
coupled field with a very particular limit in parameter space. In general, in 
dSo (and in any curved spacetime) the Klein-Gordon equation (32) includes 
the term £R(p, where R is the Ricci scalar and £ is a parameter called the 
coupling to the curvature. In this paper we have considered only minimally 
coupled fields, i.e., £ = 0, and we have compared the case £ = 0, m = with 
the case £ = 0, m — ► 0. In Ref. [35], the expectation value of the stress tensor 
in the case £ = 0, m = was compared with the same expectation value in 
the more general limit £ — > 0, m — ► 0. It was found that such a limit is not 
unique. In particular, it was shown that, if the limit is taken in such a way 
that £i?/m 2 — > — 2, the discontinuity disappears. It is reasonable to think 
that the stress tensor two-point function will exhibit a similar behavior. 

We close this subsection by mentioning that the stress tensor two-point 
function of a massless minimally coupled scalar field in dS^ was also com- 
puted in Ref. [36]. Although in qualitative agreement, the result of that 
reference differs from the results presented here because the covariant deriva- 
tives of the unit vectors were not properly taken into account there. 
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6.3 Large distances 

In the previous subsections we have studied the stress tensor fluctuations in 
the cases of small masses and strictly zero mass. In both cases we have paid 
attention to the long-distance {Z <C —1) regime. Here we investigate the 
long-distance behavior of the stress tensor fluctuations of a massive field, 
without making assumptions on the smallness of m. 

Using the properties of the hypergeometric function [33] , one sees that 
the Wightman function (40) can be expanded as 

G(n) = g+(-Z)- h + + g_(-Z)- h - + 0[Z~ ™"(h+,h-)-i^ (58) 

where 

9± (4tt)D/2 r(f-/i±) ' 1 j 

and where h± is given in (41). Now, there are two possible situations, 
depending on the size of m compared to (D — l)/2. If m < (D — l)/2, 
the parameters h± are real, with h + > h-. Then, the first term in (58) is 
subleading, and the equation simplifies to 

G(n) = g-(-Z)- h - + o(Z- min ( h +' h - + V). (60) 

After substitution in (42), one obtains the long-distance behavior of the 
stress tensor two-point function. The result is 

P,Q,T ~ Z- 2h - R~Z- 2h — x S~Z- 2h — 2 . (61) 

This is consistent with what we have seen in Sec. 6.1: if we let m — > 0, it 
follows from (41) that h- ~ m 2 /(D — 1) — ► 0, and the leading contribution 
to the stress tensor two-point function is a constant. In fact, the coefficients 
of the powers of Z that have been omitted in the equation above all vanish 
when m — > except those of S and T. Therefore, the only function that 
tends to a constant is T, as has been seen in 6.1. Besides, R decays more 
quickly than Z -1 , also consistently with that subsection. The parameter h- 
increases with the mass as long as m < {D — l)/2. Hence, the fluctuations 
decay faster with the distance as the mass increases. 

In the opposite case, namely m > (D — l)/2, the parameters h± are 
complex, and so are g±. According to Eqs. (41) and (59), we can write 

h ± = l -{D-l)±ia g±=ge ±l5 , 

where a, g and 5 are positive real numbers. In terms of these parameters, 
Eq. (58) takes the form 

G(n) = 2 5 cos[aln(-Z) - 8\(-Z)-^ + 0{Z~^). (62) 
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Substituting this equation in (42), we obtain the following result for the 
long-distance behavior of the stress tensor two-point function: 



P,Q 



,T~Z- D+1 



—D 



s ~ z 



-D-l 



(63) 



where the coefficients that have been omitted are oscillating functions of Z. 
Hence, for m > (D — l)/2, the way in which the stress tensor fluctuations 
decay with the distance does not depend on the mass of the field. This 
is different from the situation in Minkowski spacetime, where the Wight- 
man function decays with the distance like e~ m ^, and accordingly the stress 
tensor fluctuations decay like e~ 2mfJl . 

Naively, one might expect to recover the results of Minkowski spacetime 
when the Compton wavelength 1/m is much smaller than the Hubble ra- 
dius (which here is set to one). Then one would find a contradiction with 
our results above for m > (D — l)/2. However, the correct expectation is 
that one should recover the results of Minkowski spacetime when both the 
Compton wavelength and the physical distance between the two points un- 
der consideration are smaller than the Hubble radius. This is not the case 
we are dealing with here, where we assume Z <C — 1, which means that the 
physical distance is much larger than the Hubble radius. 

7 Timelike separated points 

In the previous two sections we have been discussing the properties of the 
stress tensor two-point function for spacelike separated points. Here we 
extend these results to timelike separated points. 

The Wightman function of the free scalar field in the Bunch-Davies vac- 
uum for timelike separated points, Z > 1, is [27] 



where F is the hypergeometric function, and c m ^D and h± are given in 
Eq. (41). Thus, G + has the same form as for spacelike separated points, 
except that now it depends on Z instead of Z. The biscalar Z is defined as 



where e is a positive real infinitesimal, and the minus (plus) sign is chosen 
when x is to the future (past) of x'. This =Fie term is important because the 
hypergeometric function, F(a,b;c; z), has a branch cut at z > 1 along the 
real axis. 

Introducing the biscalar fi = cos" 1 Z, the Wightman function (64) can 
also be written as 




) 



(64) 



Z(x, x') = Z(x, x') =F it-. 



(65) 



G + = G(jl), 



(66) 
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where G is the function defined in (40). From (16) and the discussion below 
that equation it follows that 



jj,(x, x') = ±fi(x, x') + e, 



(67) 



where the plus (minus) sign is chosen when x is to the future (past) of x' . 
Unlike in Eq. (16), here we are keeping the leading term that vanishes when 
e — > + , because this limit has to be taken after evaluating the hyperge- 
ometric function at Z. Recall that the geodesic distance between timelike 
separated points is imaginary, as has been emphasized under Eq. (10). 

Note that fi is not invariant under the full de Sitter group. For example, 
under the antipodal transformation it transforms as 



This is because, if x is to the future of x', then a(x) is to the past of a(x'). 
Also, due to the branch cut of the hypergeometric function, G(/2) / G(jj*) 
and, hence, the Wightman function in the Bunch-Davies vacuum for timelike 
separated points is not invariant under the full de Sitter group. This is not 
surprising: for timelike separated points, [4>(x), 4>(x')] does not necessarily 
vanish. If this commutator does not vanish, the product 4>{x)cj){x') is not a 
Hermitian operator and, in consequence, the Wightman function does not 
have to be a de Sitter-invariant biscalar. 

The stress tensor two-point function in the Bunch-Davies vacuum for 
timelike separated points is obtained by substitution of (66) in (37). Using 
(14), we find 

Fabc'd' = P(ji)n a n b n c >n d i + Q(p,)(n a n b g c i d > + n c m d >g ab ) 

+R(fi)(n a n c >g bd i + n b n d ig ac > + n a n d ig bc i + n b n c ig ad i) 



where P, Q, R, S and T are given by Eqs. (42) and (43). Therefore, the stress 
tensor two-point function for timelike separated points has the same form 
as for spacelike separated points, except that now the functions P,Q,R,S 
and T depend on fi instead of ji. 

The explicit form of these functions in the particular cases of small 
masses and massless field is obtained simply by replacing fx by fl (or, equiv- 
alently, Z by Z) in Eqs. (48)-(50) and (56). In fact, the functions of Z 
appearing in those equations do not have any branch cut. Hence, the =Fie 
term in the definition of Z is irrelevant, and those equations are valid them- 
selves also for timelike separated points. This means that, for both m = 
and m — ► 0, the stress tensor two-point function is also de Sitter-invariant 
for timelike separated points. The behavior at large separations, Z S> 1, 
is obtained similarly, by replacing Z by Z in Eqs. (61) and (63). There, 
strictly speaking we do have branch cuts (due to the non-integer exponents 



jl(a(x), a(x')) = fi*(x,x'). 



(68) 



+S{fl){gac'9bd' + gbc'Qad') + T{ji)g ab g c i d i 



(69) 



24 



in (61) and the oscillating functions that have been omitted from (63)), but 
the qualitative behavior is not affected by them. 

In order to interpret physically these results, it is convenient to express 
them in terms if the geodesic distance /j,, because there is always a geodesic 
joining two timelike separated points. For Z 3> 1, we have Z ~ exp as 
can be seen from (16). Let us concentrate on the cases m — > and m = 0. 
In the limit m — ► 0, due to the constant piece in (49), the fluctuations 
do not vanish at large separations. As we have already commented, this is 
different from what happens in Minkowski spacetime, where the fluctuations 
decay like n~ 2D . In contrast, in the massless case it follows from (57) that 
the fluctuations decay exponentially with that is, much faster than in 
Minkowski spacetime. 

8 Discussion 

In this paper we have studied the stress tensor two-point function of a quan- 
tum field in dSo- We have given general constraints on this two-point func- 
tion that hold for any quantum field theory in cLSd, assuming that the 
state of the field is de Sitter invariant and that the two points are spacelike 
separated. Next, we have derived explicit expressions for a free minimally 
coupled scalar field in the Bunch-Davies vacuum. These expressions take 
a simple form (in terms of elementary functions) in the particular cases of 
a small or vanishing mass as well as for large separations between the two 
points and general mass. (In the case of a vanishing mass the Bunch-Davies 
vacuum does not exist, but one can consider a vacuum arbitrarily close to 
it, in the sense specified in the appendix.) These simple expressions, which 
have been given in Sec. 6, are the main result of this paper. The conclusions 
of that section can be summarized as follows. When the mass m of the field 
is small, the correlations have a long range: they decay like an inverse power 
of the physical distance, with an exponent proportional to m 2 /H 2 . In par- 
ticular, in the limit m — > they do not vanish at long distances. In contrast, 
in the strictly massless case the stress tensor two-point function decays like 
the fourth power of the physical distance. There is thus a discontinuity in 
the massless limit. Such a discontinuity may seem surprising at first sight, 
but it can be easily understood as follows. In the massive case, there is a 
contribution from the mass term in the stress tensor that does not vanish in 
the limit m — ► due to the infrared divergence of the Bunch-Davies vacuum. 

By explicitly writing the dependence on the Hubble constant in our 
expressions for the stress-tensor two-point functions and taking the limit 
H — > one can recover the flat-space limit and compare with the results of 
Ref. [37] , where they were computed in the Minkowski vacuum for a general 
mass and curvature coupling. In fact, by generalizing Eqs. (32)-(43) one 
can straightforwardly (although somewhat tediously) extend our results to 
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a general value of the curvature coupling parameter. Furthermore, by an 
analytic continuation of the Hubble constant to iH and the choice of a so- 
lution of Eq. (39) with appropriate boundary conditions in that case [rather 
than the solution given by Eq. (40)] one can also obtain the stress-tensor 
two-point functions for the natural vacua in anti-de Sitter. This will give 
the general result in terms of invariant bitensors for all maximally symmet- 
ric spacetimes (Minkowski, de Sitter and anti-de Sitter). By particularizing 
to the massless conformally coupled case, one can then compare with the 
results of Ref. [38] for general conformal field theories. 

Next, we discuss how our results for the quantum correlations of the 
stress tensor can be used to obtain information on the quantum fluctua- 
tions of the gravitational field in de Sitter spacetime. Specifically, we will 
consider a perturbative quantization of the metric around a de Sitter back- 
ground within the framework of an effective field theory approach to quan- 
tum gravity [39, 40], and for simplicity we will restrict our attention to 
D = 4 spacetime dimensions. Moreover, we will neglect the contributions 
from graviton loops in order to avoid the difficulties involved in constructing 
and dealing with diffeomorphism-invariant observables [41] that arise when 
those are taken into account. This can be naturally implemented in the con- 
text of a large N expansion for N matter fields (free scalar fields in our case) 
coupled to the gravitational field, where one introduces a rescaled coupling 
constant l p = l p y/N (with l p being the usual Planck length) and expands in 
powers of 1/7V keeping l p fixed. Neglecting graviton loops amounts then to 
considering only the leading order in 1/N, i.e. the first non-vanishing order 
that contributes to any quantity of interest. For instance, the leading con- 
tribution to the two-point function of the metric perturbations is of order 
1/N. 

As a characterization of the quantum metric fluctuations at order 1/N, 
one can study the two-point correlation functions of the linearized Riemann 
tensor including the effects of matter loops. It should be emphasized that 
we are not referring to transition matrix elements but to the correlation 
functions appearing in the closed time path (also known as in-in), formal- 
ism [42, 43, 44, 45, 46]. They can be obtained either from a purely quantum 
field theoretical calculation [47] or within the framework of stochastic gravity 
[48, 49, 50, 51, 52], which provides an equivalent way of obtaining the re- 
sult at order 1/N [53, 54]. The Riemann tensor of the background de Sitter 

spacetime with appropriately raised indices is given by R ab cd = (-R/6) 5^8^ 
and its linear perturbation is, therefore, gauge invariant since its Lie deriva- 
tive with respect to an arbitrary vector vanishes [55]. In fact, both the Ricci 
tensor and the Weyl tensor (which vanishes in the background spacetime) 
are separately gauge invariant at the linear level, and together they entirely 
determine the Riemann tensor at that order. This is actually enough to 
guarantee that the connected part of their two-point correlation functions 
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at order 1/N is gauge invariant. One could alternatively consider the two- 
point function of other gauge-invariant quantities for the linearized metric 
perturbations, such as those commonly employed in the study of cosmolog- 
ical perturbations. Instead, here we consider the Riemann tensor because 
besides being a gauge-invariant object with a simple geometric meaning and 
completely characterizing the local geometry, the de Sitter invariance of its 
correlation functions can be directly analyzed, especially when using a mani- 
festly invariant formalism like the one employed in this paper. Furthermore, 
from our results for the stress-tensor two-point functions one can immedi- 
ately obtain the two-point functions of the Ricci tensor to leading order in 
the Planck length (and order 1/N), as we will see below. 

The lowest-order non-vanishing contribution to the two-point function of 
the Ricci tensor, or equivalently the Einstein tensor, is quartic in l p and (at 
order 1/N) comes from the one-loop contributions of the matter fields. Both 
from the stochastic gravity formalism [53, 52] or from the Einstein equation 
to linear order in the metric perturbations as an equation for operators in 
the Heisenberg picture (including the matter fields as well as the metric 
perturbations), one obtains the following result for the two-point function 
of the Einstein tensor at that order: 

(Gf\x) G c '/\x')) - (Gt\x)) (<«(*')) = (Svr) 2 | F%iM, 

(70) 

where G a ^ is the linearized Einstein tensor, and F a b c d , is the stress tensor 
two-point function that has been computed in this paper, with appropriately 
raised indices. Therefore, the stress tensor two-point function in dS± gives 
the leading order contribution to the two-point function of the linearized 
Einstein tensor in an expansion in powers of 1/N and l p . Note that in 
general this two-point function would receive additional contributions of 
the same order from the finite part of the local counterterms quadratic in 
the curvature which are required in order to renormalize the divergences of 
the one-loop contribution, but at this order these contributions turn out to 
vanish for the particular case of a de Sitter background. 

On the other hand, even though the two-point function of the Weyl ten- 
sor also gets contributions of order T£ and higher from matter loops, it gets 
in addition a non- vanishing tree-level contribution of order l p even in the ab- 
sence of matter loops, which was calculated in Ref. [21]. In contrast with the 
correlation functions of the Ricci tensor, the contributions of order T£ to the 
two-point function of the Weyl tensor cannot be immediately obtained from 
the stress-tensor two-point function. Nevertheless, when combined with the 
tree-level result for the correlations of the Weyl tensor in Ref. [21], our cur- 
rent result for the correlations of the Ricci tensor completely characterizes 
the quantum fluctuations of the full Riemann tensor at order 1/N and lead- 
ing order in the Planck length. It should be pointed out that the two-point 
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function corresponding to the cross-correlation of the Ricci and Weyl ten- 
sors is of order T£ (rather than l^). This can be interpreted as the absence 
of correlations between the tree-level fluctuations of the Weyl tensor, with 
an amplitude of order l p , and the fluctuations of the Ricci tensor, with an 
amplitude of order l£ . 

Note that in contrast to the existing calculations of loop corrections to 
the spectrum of cosmological perturbations in inflationary models, where the 
self-interaction of the matter fields often plays an important role [56, 57, 58], 
the only interaction vertices associated with the free fields that we have 
considered here are those corresponding to their gravitational interaction. 
Moreover, we do not consider a non- vanishing homogeneous background con- 
figuration for the scalar fields (which would give also a tree-level contribution 
to the stress-tensor two-point function [59]) and the usual treatment for the 
correlations of the scalar-type metric perturbations [60, 47, 61, 62] cannot 
be directly applied. On the other hand, the calculation for tensor metric 
perturbations should be equivalent (since the slow-roll geometry is usually 
approximated by de Sitter space in that case) but the existing results are 
either approximate [63] or restricted to the massless case [62] . Furthermore, 
considering the correlations of the Riemann tensor and expressing them 
in terms of maximally invariant tensors is very useful in order to analyze 
whether de Sitter invariance is disrupted by loop corrections. 

We close this section with a few remarks concerning possible lines of 
future research in connection with our work. An aspect that deserves further 
study is analyzing how generic our main results for light fields (namely, the 
existence of long-range stress-tensor correlations and the discontinuity of 
the massless limit) are when considering general states. For instance, for a 
wide range of masses and curvature-coupling parameters it has been shown 
that given a generic well-behaved and spatially isotropic initial state, the 
expectation value of the stress-tensor operator on a de Sitter background 
tends at late times to the same value as the expectation value for the Bunch- 
Davies vacuum [64]. Moreover, for spatially isotropic perturbations this 
conclusion remains true even when the back-reaction on the background 
geometry is taken into account [65, 66, 67, 68, 69]. This can be qualitatively 
understood as a consequence of the exponential redshift rendering negligible, 
after a sufficiently long time, the contribution to the stress tensor of any 
initial excitation of the Bunch-Davies vacuum. Hence, we expect that a 
similar conclusion will apply to the stress-tensor correlations and its main 
features, but it would be interesting to check it explicitly. Secondly, in 
future work we plan to compute the one-loop contribution (of order 1^) to 
the two-point function of the Weyl tensor, which unlike for the Ricci tensor 
cannot be immediately obtained from the stress tensor two-point function, 
and investigate the de Sitter invariance of the one-loop contributions to the 
correlation function of the full Riemann tensor and its implications. 
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A Vacuum states 

In this appendix we summarize the quantization of a free, minimally coupled 
scalar field in dSp. We pay special attention to the massless case, where 
we reexamine the issue of the non-existence of a de Sitter-invariant vacuum, 
previously addressed in Ref. [28, 35], and prove Eq. (54). For simplicity, we 
will work in D = 4 dimensions of spacetime. 

A coordinate system {77, fi} can be chosen for dS± such that the line 
element takes the form 

ds 2 = sin" 2 7][-dri 2 + dtt 2 } . (71) 

where r\ 6 (0,7r) is the so-called conformal time, $7 is a set of angles of S 3 , 
and dti 2 is the line element of S 3 . In these coordinates, the scalar field can 
be expanded in spherical harmonics on S 3 , 

cf>(r,,n) = ^[aLMULM^YLM^ + a^UlMYlMm, (72) 

L,M 

where the spherical harmonics Ylm(^) are normalized eigenfunctions of the 
Laplacian on S 3 , L = 0, ... ,00, and M represents two indices, M\ and M2, 
with Mi = —L, . . . , L and M2 = —Mi, . . . , M\. In order for the scalar 
field to satisfy the Klein-Gordon equation (32), the modes Ulm have to be 
particular solutions of the equation 

Ulm(v) ~ 2 cot 7] U LM (v) + ^1(v)Ulm(v) = 0, (73) 

where 0^(77) = L(L + 2) + m 2 sin -2 ?7. If these particular solutions satisfy 
the Wronskian condition 

Ulm{v)U1m(v) ~ UI m {v)Ulm{v) = * sin 2 77, (74) 

then the canonical commutation relations between the field and its conju- 
gate momentum imply that a} LM and aiM are creation and annihilation 
operators, 

[ a LM,a^ L ' M i] = Sll'&mm'- (75) 
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The vacuum |0) is then defined as the state annihilated by all the annihila- 
tion operators, 

a LM |0)=0 VL,M. (76) 

In spite of the Wronskian condition (74), there is freedom in choosing the 
modes Ulm- Different choices give rise to different vacua. The Bunch-Davies 
vacuum, \0)bd, is associated with the following choice [70, 27] 

UEm(v) = ALsm^niPX-cosr!) - -Q*(- cost?)], (77) 

where P^ and are Legendre functions on the cut, A = (9/4 — m 2 ) 1 / 2 , 
v = L + 1/2, and the normalization constants are 



r(L- A + 3/2) 



2 [r(L + A + 3/2) 



1/2 

(78) 



The Bunch-Davies vacuum is de Sitter-invariant. Note that this state is 
only defined for massive fields, m > 0, because Aq does not exist for m = 0. 
In fact, in the massless case there is no de Sitter-invariant vacuum, as was 
shown in Ref. [28]. One can instead consider the vacua \0} a , associated with 
the following modes [34] : 

US M (v)= UmUEgirj) forL>0 

[/ a (r ? ) = -m(r ? -^sin2r ? -|) + ^, (79) 

where a £ (0,oo). These vacua are invariant under the 0(4) subgroup of 
0(4, 1), i.e., the subgroup of isometries that leave invariant the hypersurfaces 
of constant ry. They are also invariant under time reversal. (These vacua 
should not be confused with the so-called alpha vacua, which are defined for 
massive fields and whose Wightman function diverges for antipodal points.) 

In what follows we concentrate on the massless case. As was pointed out 
in Ref. [35], the non-existence of the Bunch-Davies vacuum when m = is 
better understood by looking at the example of the harmonic oscillator in 
Quantum Mechanics. The general solution of the equation of motion of a 
harmonic oscillator of frequency to can be written in the form 

q(t) = U(t)a + U*(t)a\ (80) 

where and a are creation and annihilation operators, and U(t) is a par- 
ticular solution of the equation of motion. If it is chosen as 

U(t) = -ge"**, (81) 



then the ground state \ip) of the harmonic oscillator satisfies a\ip) = 0. Now, 
U{t) does not exist when w = 0. This is equivalent to what happens to 
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Uq when m = 0. But in this case the interpretation is clear: a harmonic 
oscillator of zero frequency is a free particle. And the ground state of the free 
particle, which should be annihilated by the momentum operator, p\ip) = 0, 
does not exist, because a constant wave function cannot be normalized. 
Similarly [35], the wave functional of the Bunch-Davies vacuum becomes 
non-normalizable in the limit m — > 0. 

Even though the ground state of the free particle does not exist in Quan- 
tum Mechanics, there are states which are "arbitrarily close" to it. An ex- 
ample is a Gaussian wave function with very large width. Similarly, one can 
expect to find states for the massless field in dS± which are "arbitrarily close" 
to being de Sitter-invariant. This is indeed true, in the following (restricted) 
sense: even though the vacua |0) a are not de Sitter-invariant, there is a cer- 
tain class of observables whose expectation values are de Sitter-invariant in 
the limit a — > 0. 

Let us see this. For any mass m, and in any vacuum state |0), the 
Wightman function of the free scalar field in dS± is expanded in terms of 
the modes as 

G + (r ] ,n; V ',n') = Y,YLM(n)YZM(rt)U LM (r))U* LM (r l > ). (82) 

L,M 

If we want to compute the Wightman function, G+, of a massless field in 
the state \0) a , we have to replace Ulm by U" M in this equation. If, instead, 
we want to compute the Wightman function, G~^ D , of a massive field in the 
Bunch-Davies vacuum, we will use U^m- Let us consider the limit a — > 
for a massless field in the state |0) a . It follows from (79) that the zero 
mode Uq diverges in this limit, and so does the Wightman function G+. 
However, the derivative of the zero mode, Uq, vanishes. Therefore, d^d^G^ 
is finite in the limit a — > and, besides, the term with L = does not 
contribute to it. On the other hand, consider the limit m — > for a massive 
field in the Bunch-Davies vacuum. It follows from (78) that the zero mode 
Uq D diverges. This is the origin of the infrared divergence of G~g D found in 
Sec. 6.1. However, the derivative of the zero mode, Uq D , vanishes. Hence, 
drjdjjiG^p is finite in the limit m — > 0, and the term with L = does not 
contribute to it. Using Eqs. (79) and (82) we find 

lim d ri d n iGt = lim d ri d n iG% n . (83) 

Finally, the spherical harmonic loo is independent of S7. Therefore, whenever 
we take derivatives of the Wightman function (82) with respect to £1 or f2' 
the term with L = will not contribute, and a result similar to the above 
equation will apply. From all this we conclude that 

lim V a V 6 ,G+ = lim V a V 6 ,G+ D . (84) 
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Since the right-hand side is de Sitter-invariant, this equation implies that 
V a Vfe'Gj is de Sitter-invariant in the limit a — ► 0. The expectation value 
in the vacuum |0) Q of any observable that only depends on derivatives of 
the field can be expressed in terms of V a V;/G+ , because |0) Q is a Gaussian 
state. In consequence, such expectation value will be de Sitter-invariant in 
the limit a — > 0. Hence, we conclude that, even though the states |0) a are 
not de Sitter- invariant, the expectation values of certain observables can be 
taken arbitrarily close to a de Sitter-invariant value by choosing a sufficiently 
small. 

Note that the conclusions above apply only to a certain class of observ- 
ables, not all of them. For example, V a G+ does not satisfy an equation 
like (84). The reason is that, when computing this bitensor using Eq. (82), 
the contribution from the term with L = is finite, but it does not vanish. 
Hence, one cannot use the argument that led to Eq. (84). It is also worth 
noting that |0) Q is not the only family of states satisfying (84). For instance, 
consider the family of coherent states \z) a , defined by a,LM\z)a = for L > 
and ao\z) a = z\z) a , where z is a complex number, and the annihilation op- 
erators are those associated with the 0(4) invariant vacua |0) Q . It is not 
hard to see that this family of states also satisfies Eq. (84) for all values of 
z. 
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